BCS-BEC crossover on the two-dimensional honeycomb lattice 
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The attractive Hubbard model on the honeycomb lattice exhibits, at half-filling, a quantum critical 
point (QCP) between a semimetal with massless Dirac fermions and an s-wave superconductor (SC). 
We study the BCS-BEC crossover in this model away from half-filling at zero temperature and show 
that the appropriately defined crossover line (in the interaction-density plane) passes through the 
QCP at half-filling. For a range of densities around half-filling, the "underlying Fermi surface" of 
the SC, defined as the momentum space locus of minimum energy quasiparticle excitations, encloses 
an area which changes nonmonotonically with interaction. We also study fluctuations in the SC and 
the semimetal, and show the emergence of an undamped Leggett mode deep in the SC. Finally, we 
consider possible implications for ultracold atoms in optical lattices and the high temperature SCs. 

PACS numbers: 03.75.Kk, 03.75.Ss, 74.20.-z 



The close connection between the Bardeen-Cooper- 
SchriefFer (BCS) theory of superconductivity and the 
phenomenon of Bose-Einstein condensation (BEC) is 
clearly revealed by studies of the BCS-BEC crossover 
problem Q . One can experimentally access this crossover 
in harmonically trapped atomic gases by tuning the s- 
wave atom-atom scattering length 0= using a magnetic 
field induced Feshbach resonance 0j Hi S IE IE • The 
low temperature phase of these atomic gases crosses over 
from a BCS state for Os < to a molecular BEC for 
as > 0. Experiments H and theory [IE El 113 
have also begun to address the issue of strongly interact- 
ing fermionic atoms in optical lattices. The ETH mea- 
surements IE) of the population distribution in different 
bands upon tuning through a Feshbach resonance were 
explained partially by focussing on a single well in the 
lattice EE- I'^ addition, it was conjectured that this sys- 
tem could show a band insulator to SC transition with 
increasing interaction (IQ] at low temperature; however 
multiband effects EE complicate a microscopic analysis 
of this problem. 

Motivated by exploring a simpler model which displays 
such a transition and can be realized using cold atoms, 
we study the attractive Hubbard model, 
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on a two-dimensional (2D) honeycomb lattice, with hop- 
ping amplitude t.y = t to nearest-neigbor sites and 
tij = t' to next-nearest neighbors (see Fig.l). Nonin- 
teracting fermions {U = 0) at half-filling (one fermion 
per site on average) on this lattice form a semimetal, 
with the "Fermi surface" shrunk to Fermi points. This 
semimetal behaves like a band insulator in some respects, 
displaying a vanishing density of states at zero energy and 
two bands (a "conduction" and a "valence" band) which 
touch at two inequivalent Fermi points, ±K, in momen- 
tum space as shown in Fig.l (the two bands arise from 
having two sites per unit cell on the honeycomb lattice). 
We show that turning on interactions in this model leads 
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FIG. 1: Left: Honeycomb lattice showing basis vectors and the 
two sites per unit cell, as well as the hopping matrix elements 
in model (1). Right: Quasiparticle dispersion with Dirac 
nodes in first Brillouin zone of the noninteracting semimetal. 



to rich physics associated with a semimetal-SC transition 
at half-filling and a BCS-BEC crossover away from half- 
filling when U becomes comparable to the kinetic energy 
^ t. Bands other than those in model can, however, 
be ignored so long as they are separated by interband 
gaps > [/, i. 

Compared to earlier studies of the BCS-BEC crossover 
EE I our following results are unusual, (i) The conven- 
tional BCS-BEC crossover is unrelated to any quantum 
phase transition (QPT), describing the smooth evolution 
from a weakly paired to a strongly paired SC state. By 
contrast, model (1) at half-filling exhibits a QPT be- 
tween a semimetal and a SC at a critical interaction, 
U = Uc- (u) Away from half-filling, it displays a BCS- 
BEC crossover. As shown in Fig. 2 (a), the appropri- 
ately defined BCS-BEC crossover line in the interaction- 
density plane passes through the quantum critical point 
(QCP) at half-filling, suggesting that the finite temper- 
ature normal state near the crossover could be strongly 
infiuenced by this QCP, especially for small density devi- 
ations from half-filling ("small doping"), (iii) For small 
dopings, the "underlying Fermi surface" (UFS) of the SC 
|lJ> ll^ ) defined as the momentum space locus of mini- 
mum energy quasiparticle excitations, exhibits nonmono- 
tonic behavior with increasing U in contrast to its mono- 
tonically shrinking area in the continuum EE- (i^) The 
fluctuations in the SC phase of model (1) are like those 
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of a two-band SC [ig and support an undamped Leggett 
mode deep in the BEC regime (Fig. 3(a)). (v) In con- 
trast to a Fermi liquid, long wavelength SC fluctuations 
in the semimetal have a damping proportional to energy 
(Fig. 3(b)) due to the Dirac fermions. We discuss the rel- 
evance of some of these results the high temperature SCs 
in the concluding section. 

Model and mean-field theory: We write the partition 
function of model Q as Z = J Viptl; exp{— S) , with 

- Uj2t,a'>Pa,h^i'^)'^a,t,liT)'>Pa,r,liT)->pa,r,\{T)^ where i 

refers to a site on the triangular lattice, and a — 
1,2 labels the two sites within the unit cell (sublat- 
tice index). The matrix elements /ia/3(k) are given by 
/iii(k) = /i22(k) = Xk, and /i2i(k) = K^Ck) = 7k 
with Xk = — 2t'(cosk • a + cosk • b + cosk • c) — /i and 
7k — ~t{l + e*'' '' + e^'*' °). Here a,b are unit vectors 
shown in Fig.l, and c ^ d + b. We introduce [l^ 
Hubbard-Statonovich fields Aa^i{T) to decouple the in- 
teraction term in the pair channel. The mean-field theory 
(MFT) of the SC is obtained by setting A„,j(r) = Aq, 
integrating out the fermions, and extremizing the action 
with respect to Aq. This leads to the gap equation 

Here ly labels the two bands ^ \/^±(^) + A§ 

is the excitation energy of Bogoliubov quasiparticles in 
MFT with ^±(k) = Xk ± |7k|, and N is the number of 
sites on the honeycomb lattice. Using dJ-Mp/dfi^—Ne, 
the fermion density n = N^/N is 

^=1-]^ -^tanh(^). (3) 

Figs.2(a,b) show the MF values of Aq, ^ obtained by nu- 
merically solving Q and (jSJ over a range of densities and 
interaction strengths. 

At half-filling {n = 1) and for t' < t/3, the noninteract- 
ing {U = 0) ground state is a semimetal with Aq = and 
(x = 3i'. The single-particle spectrum consists of mass- 
less Dirac fermions, centered around K = (±47r/3, 0) and 
dispersing at low energy as ^±(k) ~ ±?;i?|k — K|, with 
a Fermi velocity vp — iV3/2. This leads to a linearly 
vanishing density of states (DOS) of single-particle exci- 
tations at low energy, N{uj) = -^tUj per spin, which ren- 
ders the semimetal stable to weak attractive interactions 
in contrast to a Fermi liquid. At large enough interac- 
tions U > Uc, the semimetal becomes unstable to an or- 
dered state, similar to that in a model studied by Nozieres 
and Pistolesi For i' = 0, both SC and charge den- 
sity wave solutions are degenerate due to a special SU(2) 
symmetry. This degeneracy is lifted at nonzero t' , and 




FIG. 2: (a) (Color online) Plot of the mean field gap, Ao/f, 
versus n and U/t (for t' = —0.15t). White horizontal line at 
n = 1 is the semimetal phase terminating in a semimetal-SC 
transition at Uc/t « 2.13. Dashed lines indicate the BCS- 
BEC crossover discussed in the text. Thick black line indi- 
cates the parameters at which the "underlying Fermi surface" 
(UFS), discussed in text and panels (b) and (c), shrinks to a 
point. The light arrow indicates a possible trajectory of the 
Hamiltonian discussed in the text, (b) Mean field chemical 
potential ^/t versus U/t for indicated fillings {t' = — 0.15t). 
Dashed line at /i = 3t' is a guide to the eye. (c) Schematic 
"underlying Fermi surface" (see text), roughly circular loci 
around the nodes ±K, drawn for parameters indicated by 
arrows in panel (b). 

we find that the SC state has lower energy justifying our 
Hubbard-Stratonovitch decoupling in the pairing chan- 
nel. The mean field gap Aq ~ (C/ — Uc) for U close to Uc 
while for large U, Aq ~ U/2. We fix t' = -0.15i, where 
Uc ~ 2.13t, and focus here only on the SC order and its 
fiuctuations. (For t' — 0, MF, variational, and Monte 
Carlo numerics yield estimates Uc ^ 2t — 5t l2l|.) 

Away from half-filling, with a doping 6n = n — I, the 
noninteracting ground state has small Fermi pockets cen- 
tered around ±K, and is unstable to SC for an arbi- 
trarily small U. The DOS scales as ^/\^n\/vF leading 
to Aq ~ \/\Sn\ expl—gvp/iU \/\Sn\)] for U ^ VF\/\Sn\ 
(with g'^ = 167r/\/3), where VF^/\6n\ acts like an effec- 
tive Fermi energy. This weak coupling SC evolves into a 
BEC regime at larger U. 

(i) BCS-BEC crossover: A mean field estimate of the 
crossover to the BEC regime is obtained by requiring 
that the SC gap Ao{U,6n) be larger than either (i) the 
Fermi energy of the doped carriers in the noninteracting 
limit /i(C/ = 0, 6n) — 3t' or (ii) the chemical potential 
tiiU^n) - 3t', where 3t' = n{U = 0,Sn = 0). Fig.2(a) 
shows this crossover line plotted using criterion (i); how- 
ever both definitions yield the qualitatively similar re- 
sult that the BCS-BEC crossover line (defined away from 
half- filling) passes through the QCP (at half-filling). If 
the Hamiltonian describing the doped system lies on a 
trajectory such as shown in Fig. 2(a), with a density de- 
pendent U, then the system will undergo a BEC-BCS 
crossover for parameters which lie close to the QCP; the 
finite temperature normal state near the crossover is then 
expected to be strongly infiuenced by this QCP. 
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(ii) Underlying Fermi surface: We define the "underlying 
Fermi surface" (UFS) of the SC as the locus of points 
around ±K where the quasiparticle excitation gap is a 
minimum [l^ . At weak coupling, this UFS is identical 
to the FS of the doped metallic phase, which is approx- 
imately a circle centered around ±K. Turning on inter- 
actions at small doping leads to a strong change in /i. In 
fact, for a range of doping Sn < (for t' < 0), fi goes 
from the "valence" into the "conduction" band leading 
to the evolution of the UFS as shown in Fig.2(c). Even- 
tually, at a large enough gap, the UFS begins to decrease 
[2^ as in the continuum case The UFS is thus not 
related in any way to the Fermi surface of the underlying 
noninteracting metal |l4l For example, at a certain 
interaction value for low doping we find ^ = 3t'; since the 
chemical potential then sits at the node of the Dirac dis- 
persion it leads to a point-like UFS as seen in Fig. 2(c). 
For this value of U, the strongly interacting doped SC 
would appear, if naively viewed as a weak coupling SC, 
to have descended from an undoped semimetal. 
Superfluid density, Estimate of Tc-' Going beyond MFT, 
we expand Acq(T) = Aq -I- AQq(T) around its mean 
field value and integrate out the fermions to arrive at 
the effective action for order parameter fluctuations, 
•^RPA = J2q^l^q^q^ to second order in Aa,q. Here 
q = (q,iw„), and A, = (Ai,^, A| Aj,,, A^ )T de- 
fines the vector of order parameter fluctuations on the 1, 2 
sublattices. We further separate the fluctuation into its 
amplitude (real) and phase (imaginary) components 
Aa^q = Ao[ra,q + i(pa,q], and integrate out the (gapped) 
amplitude fluctuations to obtain the effective 2x2 matrix 
action for phase fluctuations 



iV*l,q, V*2,q) 



^q "q 
'*q Uq 



(4) 



Diagonalizing this leads to two eigenmodes, the Gold- 
stone mode ipG and the Leggett mode ip^, with S[ip\ ~ 

(/3/2) [K - \Vq\)WG,q? + {Uq + \Vq\)WL^q\'] . Thc 

zero temperature superfluid density, DsijS), is obtained 
by examining the Goldstone mode action in the limit 
bj — Q. In this static limit, amplitude and phase fluctua- 



tions decouple, and u„ 



2A| 

U 



2N ^ 



Ua,Va 



COS 



Vq simplify to 



2N ^ 



COS 



k-q 
k-q^ 7k7k-q 



^k 



k-q 



l7k||7k- 



(5) 
(6) 



where z/, a 



±, and 6*1^ are defined by sin20{^ = AjE^^ 
— ^^(k)/£'j^. We numerically extract the 
zero temperature superfluid stiffness, -0^(0), by identify- 
ing - It^ql = (\/3i?,(0)/2)q2 for |q| 



and cos20j^ 



0. 



The SC-normal transition temperature at weak cou- 
pling is mainly determined by the vanishing of Aq due to 



thermally excited quasiparticles. The mean field transi- 
tion temperature T^, where Aq 0, is obtained by gener- 
alizing and solving the mean field gap and number equa- 
tions at nonzero temperature [23|. At strong coupling, 
Tc is governed by the superfluid stiffness. We estimate 
the Kosterlitz-Thouless (KT) transition temperature as 
Tjjrp = 7rZ?s(0)/2. The SC transition temperature is then 
approximately given by Tc = min(T°, TjJx)- Fig-3(a) dis- 
plays the crossover region between the BCS (Tc ~ r°) 
and BEG (Tc ^ T^t) regimes. We also find that the 
highest transition temperature T^'^^ ^ O.li. 
Leggett mode, Goldstone mode: In order to study collec- 
tive modes in the SC, we retain time-dependent fluctu- 
ations around the MFT, which couples amplitude and 
phase fluctuations leading to a 4 x 4 fluctuation ma- 
trix. Examination of the four eigenmodes reveals that 
one mode is gapless and corresponds to a linearly dis- 
persing Goldstone mode. The remaining three modes 
are gapped: two of these correspond to amplitude fluc- 
tuations while the third mode corresponds to a Leggett 
mode ^3 (relative fluctuations of i/j within a unit cell) 
of this "two-band" superfluid 0|. Integrating out the 
amplitude modes, leads to complicated expressions for 
Uq,Vq which will be presented in detail elsewhere [2^ . 
Using these, we have numerically evaluated the Gold- 
stone mode velocity, cq [i^l, as well as the gap to the 
Leggett mode, wl, as a function of U and n. At weak 
coupling, the Leggett mode at q = has weight only at 
energies above 2Ao and is strongly damped by the two- 
quasiparticle continuum in the SC. With increasing U , 
however, it emerges undamped (below the continuum) 
in the regime indicated in Fig. 3(a). a;L/2Ao decreases 
monotonically with increasing U in this regime. The 
Leggett mode can be excited by making small oscillations 
in opposite directions of the optical lattice potential on 
the two sublattices ;22]. 

Fluctuations in the semi-metal: The Gaussian action in 
the semimetal can be derived by studying fluctuations 
around Ao = 0. We again separate the fluctuations into 
real and imaginary parts (these no longer have the mean- 
ing of 'amplitude' and 'phase' fluctuations) and obtain 
the propagator (r* ^r^.q) — (<y5* q^a,q) from the fluc- 
tuation matrix 22]. In Fig.3(b), the spectral function 
A(q, u) = Im(r* ,jrc^g) is plotted for various values of q 
aiU — 2.1t. The lineshape of A(q, w) can be understood 
by analyzing the DOS of the decay channel for SC fluc- 
tuations in the semimetal, viz. two-particle excitations 
with total momentum q. At q = 0, this is proportional 
to the single particle DOS in the semimetal and increases 
linearly with energy. For small but nonzero q there is a 
threshold vp\q\ above which damping onsets. For large 
momenta, the threshold again vanishes at the wavevector 
connecting the Dirac nodes. Thus, the long wavelength 
SC fluctuations, which have a gap Eg that vanishes as 
U ^ Uc, are critically damped (i.e., have damping ~ 
energy) if |q| < Eg/vp. For larger q the mode begins 
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viability of realizing this model with ultracold atoms in 
an optical lattice, it would be interesting to explore its 
phase diagram in experiments. 
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FIG. 3: (Color online) (a) Plot of the SC transition tem- 
perature, Tc = min(r°,rKT) (see text). Shaded region is 
the crossover regime between Tc ~ T° ("BCS") and ~ 
^KT ( "BEC" ) . Region marked "Leggett" , beyond the dashed 
line, indicates parameters where an undamped Leggett mode 
emerges. Solid horizontal line at n = 1 is the semimetal termi- 
nating at a QCP. (b) Spectral function (in arbitrary units) of 
SC fluctuations in the semimetal at i7 = 2.1t. Small q fluctu- 
ations appear critically damped, with the damping threshold 
being strongly q-dependent. Large q fluctuations can propa- 
gate as undamped modes when close to the QCP. 

to emerge from the two-particle continuum and sharpens 
(see Fig. 3(b)). For large enough q it can propagate as an 
undamped mode for U sufficiently close to Uc- 
Discussion: The honeycomb lattice attractive Hubbard 
model exhibits a quantum phase transition between a 
semimetal and a SC, and a BCS-BEC crossover at 
low doping which lies near the transition. This model 
could be studied experimentally using two-component 
fermionic atoms in a honeycomb lattice potential which 
can be generated using three blue-detuned laser beams 
[23! . Our work is also of broad relevance to the high tem- 
perature cuprate SCs. The cuprate SCs have a d-wave 
order parameter and descend from a Mott insulating anti- 
ferromagnet Nevertheless, studies of the BCS-BEC 
crossover problem in lattice s-wave SCs P, 0| clarified 
the physics of SCs with a large pairing gap and a small 
superfluid density, and were relevant to the understand- 
ing aspects of the pseudogap regime. There are similar 
connections between our model and the cuprates. (a) 
Superconductivity in the cuprates is widely viewed as re- 
sulting from doping a spin liquid insulator |24| — a model 
relevant to the undoped spin liquid is the staggered flux 
phase jijl which, in mean field theory, has nodal Dirac 
fermions similar to our undoped semimetal. (b) The SC 
crosses over with increasing doping from a strongly corre- 
lated (BEC-like) to a weakly correlated (BCS-like) state 
0, with the normal state above Tc at optimal doping 
displaying proximity to an as-yet-unknown QCP. (c) Re- 
cent photoemission experiments in underdoped cuprates 
appear to find an anomalous metallic phase with a point- 
like FS as the normal state of underdoped cuprates [2^. 
Our results on the crossovers in a SC state close to a 
semimetal-SC QCP, and the evolution of the UFS in the 
crossover regime could shed light on observations (b),(c) 
in the cuprates. Given these connections, and also the 
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